Long-Time Tails, Weak Localization, and Classical and Quantum Critical Behavior 
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An overview is given of the long-time and long-distance behavior of correlation functions in both 
classical and quantum statistical mechanics. After a simple derivation of the classical long-time 
tails in equilibrium time correlation functions, we discuss analogous long-distance phenomena in 
nonequilibrium classical systems. The paper then draws analogies between these phenomena and 
similar effects in quantum statistical mechanics, with emphasis on the soft modes that underly long- 
time tails and related phenomena. We also elucidate the interplay between critical phenomena and 
long-time tails, using the classical liquid-gas critical point and the quantum ferromagnetic transition 
as examples. 
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I. INTRODUCTION 

It has been known for some time that equilibrium time 
correlation functions gencrically show a power-law tem- 
poral decay in the limit of long times. If this behavior is 
not obviously related to any conservation law, or to soft 
modes in the system, long-time tails (LTTs) is the term 
traditionally applied to this phenomenon. Boltzmann- 
type kinetic equations usually predict an exponential de- 
cay in time, which is why the discovery of LTTs came 
as a considerable surprise. An example is the velocity 
auto-correlation function in a hard-sphere fluid, which is 
shown in Fig. ^ In accord with the Boltzmann-Enskog 
equation, initially the decay is exponential, with a time 
scale given by the mean-free time to between particle 
collisions. However, for longer times, t >> to, it de- 
cays algebraically, as f"'^/^, and more generally as i""*/^ 
in d-dimensions. This striking result was first observed 
in computer simulations by Alder and Wainwright,EI and 
then understood theoretically by Dorfman and Cohen,El 
and by Ernst, Hauge, and van Leeuwenfl The velocity au- 
tocorrelation function is of particular importance because 
it determines the frequency-dependent self-diffusion co- 
efficient D{oj) via 



D(uj) 



dt e" 



(v(t).v(O)) 



(1.1) 



Here v(i) is the velocity of a tagged particle at time t, 
and the angular brackets denote an equilibrium ensemble 
average. The slow algebraic decay of the velocity auto- 
correlation function implies a low-frequency behavior 
D{u: ^ Q)/D{id) = 1 - (iLjto)(^-2)/2 ^ ^ ^ ^ 



D{ijj 0)/D{0) = 1 - ln{iLuto) + 



(2 < d < 4) , 
(1.2a) 

. (d=2) , 

(1.2b) 
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FIG. 1. Normalized velocity autocorrelation function 
poit) = {v(t)-v(0))/{v^(0)) as a function of the dimensionless 
time t* = t/to, where to is the mean-free time. The crosses 
indidate computer results obtained by Wood and Erpenbeck, 
Ref. ^ for a system of 4000 hard spheres at a reduced density 
corresponding to 1//Vb = 3, where V is the actual volume and 
Vb is the close-packing volume. The dashed curve represents 
the theoretical curve pn 

(t) = an (t*)-^/^ The solid curve 
represents a more complete evaluation of the mode-coupling 
formula with contributions from all possible hydrodynamic 
modes and with finite-size corrections included, Ref. a. From 
Ref. |. 



with > 0. This result indicates that in d < 2 the 
zero- frequency diffusion coefficient does not exist, which 
in turn implies that the ordinary local hydrodynamic de- 
scription of a fluid is not valid in low dimensions. This 
conclusion is indeed correct, as detailed investigations 
show.l3 Note that in this classical case transport in low 
dimensions is faster than diffusive. More generally, it is 
known that all the transport coefficients in a fluid are 
given by time integrals over appropriate time correlation 
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functions, and thai, all these time correlation functions 
have similar LTTsB 

In contrast to this behavior of time correlation func- 
tions, the spatial decay of static correlations at large dis- 
tances in classical systems in equilibrium is generically 
indeed exponential. An exception is the behavior at crit- 
ical points. It was known long before the discovery of 
LTTs that at these isolated points in the phase diagram 
correlations decay algebraically in both time and spiacje, 
i.e., like power-laws with universal critical exponents.LTtil 
Algebraic decay implies scale invariance, i.e., the corre- 
lation functions are generalized homogeneous functions 
of space and time. Correlations that exhibit scale invari- 
ance in an entire region of the phase diagram, rather than 
just at isolated paiHts, are said to exhibit generic scale 
invariance (GSI).eH13 In this language, the LTTs are an 
example of GSI in the time domain. 

A natural question is what happens to the correla- 
tion functions that exhibit GSI as a critical point is ap- 
proached, and how these correlations influence the criti- 
cal behavior itself. We will see below that the mechanism 
that causes the LTTs in the entire fluid phase typically 
becomes amplified near critical points and leads to non- 
trivial (i.e., non-van Hove) critical dynamics, or critical 
singularities in various transport coefficients. Violation 
of the van Hove picture for the slowing down of critical 
fluctuations became apparent when the thermal conduc- 
tivity of fluids was found to diverge at the critical point, 
while the viscosity of liquid mixtures exhibited a weak 
divergence at the consolutc point These singularities 
were later understood to be due to mode-coupling e£= 
fects by Fixman,0 and then by Kadanoff and Swifttj 
and by Kawasaki. tj Here mode-coupling means the cou- 
pling of products of slow hydrodynamic modes, or excita- 
tions, that are in some sense orthogonal to single hydro- 
dynamics modes, to the single modes whose correlations 
form the correlation functions in question. Historically, 
the LTTs were discovered after these critical singulari- 
ties were observed, and it turneci|Out that they, too, are 
caused by mode-coupling effects.H 

In a classical system in a nonequilibrium steady state 
(NESS) the situation is still more complicated. As we will 
discuss at the end of the first part of the present paper, in 
this case the same mode-coupling correlations that lead 
to the equilibrium LTTs and to the critical singularities, 
also lead to GSI in spatial correlations. 

In recent years there has been a separate but analogous 
development rin the study of electronic systems at zero 
temperatureJl3 In this case the mode-coupling effects are 
stronger than in classical systems, leading to even more 
dramatic effects. For example, they lead to the local- 
ization of electrons in disordered two-dimensional sys- 
tems, even for arbitrarily weak disorderJla In addition, 
because of the coupling of the electron dynamics to the 
various electronic correlation functions, and because of 
the coupling between statics and dynaipics that is inher- 
ent in quantum statistical mechanicsji3 it turns out that 
in zero-temperature systems, GSI exists in both spatial 



and temporal correlations even in equilibrium. In disor- 
dered systems these many-body quaatum phenomena are 
known as weak-localization effects,c3 since they can be 
considered precursors to aayUiderson or Anderson-Mott 
metal-insulator transition,EirE2l but analogous, if slightly 
weaker, effects also occur in clean quantum systems. In 
the second part of this paper we explain these phenom- 
ena, making connections with the GSI in classical sys- 
tems. We then explore how these generic correlations in- 
fluence, and are themselves modified by, critical behavior 
at quantum phase transitions. The latter are defined as 
phase transitions that occur at zero temperature and are 
triggered by some non-titewHal control parameter, like 
pressure or composition BiiEa The metal-insulator tran- 
sitions mentioned above are examples of quantum phase 
transitions. 

As noted above, in the classical case the basic source 
of the power-law decay characteristic of GSI is two-fold. 
First, soft or massless modes must exist. Their softness 
can arise either from some basic underlying conservation 
law, or they can be Goldstone modes, i.e., arise from a 
spontaneously broken continuous symmetry. Second, the 
theory must in some sense be nonlinear, so that these 
soft modes couple to the physical observables. It is this 
coupling that is missing in simple theories, and this omis- 
sion results in the erroneous prediction of exponential 
decay. In both classical and quantum systems such cou- 
pling mechanisms exist, and in this way the GSI effects 
found in both systems are related. 

The above discussion illustrates that it is incorrect to 
view LTT effects, as it is sometimes done, as only lead- 
ing to small corrections to transport coefficients, with no 
deeper importance. To the contrary, as we have pointed 
out, they are responsible for a myriad of effects. Not only 
do LTTs cause the leading frequency and, in the quan- 
tum case, temperature corrections to the transport coef- 
ficients, but there are numerous cases where their effects 
are greatly magnified, and they dominate the physics 
both qualitatively and quantitatively. Examples include, 
systems in reduced dimensionalities, systems near critical 
points, systems in Goldstone phases, and nonequilibrium 
systems. One of the goals of this paper is to illustrate how 
these effects come about, and how they are connected. 

We will proceed as follows. In Sec. || we present the 
dynamical equations that describe fluctuations in a clas- 
sical fluid. We use these equations to first describe the 
LTTs in a classical fluid, and then to describe critical 
dynamics near a liquid-gas critical point. We conclude 
this section by describing the GSI that occurs in spatial 



correlations in a nonequilibrium fluid. In Sec. Ill we 
consider the weak-localization effects in disordered elec- 
tron systems, which are quantum analogs of the classical 
LTTs. Because in quantum statistical mechanics these 
dynamical correlations couple to the static correlations, 
they lead to GSI in spatial correlations even in equilib- 
rium systems, as long as they are at zero temperature. 
We then discuss how this GSI is important for describ- 
ing quantum phase transitions such as the ferromagnetic 
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transition in itinerant electron systems at zero tempera- 
ture. This amounts to a discussion of how the quantum 
LTTs affect quantum phase transitions. We also discuss 
how nonequilibrium effects lead to even longer-range cor- 
relations in quantum systems at zero temperature. We 
conclude with a discussion of our results in Sec. 



II. GSI AND PHASE TRANSITIONS IN 
CLASSICAL SYSTEMS 

A. Fluctuating hydrodynamics 

To be specific, let us consider a classical fluid as an ex- 
ample of a finite temperature, or classical, statistical me- 
chanics system. In general we are interested in the long- 
time and large-distance behavior of systems, so we will 
always be concerned with the identification of the rele- 
vant soft, or hydrodynamic, variables that determine this 
behavior. In our case they are the conserved variables of 
mass density p, momentum density g = pu, with u the 
fluid velocity, and energy density e, and the relevant low- 
frequency and long-wavelength dynamical equations are 
the Navier-Stokes equations. To describe fluctuations, a 
Langevin force is added so that the equations areE3 



(2.1a) 



+ (C--^^^^)'5a/3(V-u)+P„^] , (2.1b) 

pT(9t-|-u- V)s = V„(AVc.T + fc) . (2.1c) 

Here p is the pressure, T the temperature, rj the shear 
viscosity, C, the bulk viscosity, and A the thermal con- 
ductivity, and a summation over repeated indices is im- 
plied. We have chosen to use the entropy density s, rather 
than the energy e, as our hydrodynamic variable, and in 
Eq. ( ^.Ic ) we have neglected a viscous dissipation term 
that represents entropy production since it is irrelevant 
to both the leading LTTs, and to the leading critical dy- 
namics singularities. The Langevin forces Pap and are 
uncorrelated with the initial velocity, and satisfy 

(PQ/3(x,t)Ppi.(x',i')> ^"^kBT lliSap^Sp^ + SauSp^j) 

(2.2a) 

(fc(x,i)g/3(x',i')> =2fcBAT2,5,^J(x-x')<5(t-t') , 

(2.2b) 

(P„/3(x,0<ZM(x',t')> =0 , (2.2c) 

where fee is Boltzmann's constant. The above equations 
can be derived in a number of ways and are known to 
exactly describe the Ipag-wavelength and low-frequency 
fluctuations in a fluidHD 



B. Long-time tails 

To illustrate the LTTs we choose a slightly different ex- 
ample than the self-diffusion coefficient discussed in the 
Introduction, namely, the shear viscosity of a classical 
fluid. We present a simplified calculation of this trans- 
port coefficient, and then discuss more general results. 
For simplicity, we assume an incompressible fluid. The 
mass conservation law is then trivial. 



V • u(x,i) = 



(2.3) 



and the momentum-conservation law is a closed PDE for 



dtUa -h (U • V)Ma = -^a{p/p) + V^'^U^ + "^jdPap/p ■ 

(2.4) 

Here v — ry/p is the kinematic viscosity, which we assume 
to be constant. The only role of the pressure term is to 
ensure that the incompressibility condition is satisfle d. I n 
fact, it can be eliminated by taking the curl of Eq. (|]J), 
which turns it into an equation for the transverse velocity, 
see Ref. ^ and below. The cause of the LTTs is the 
coupling of s low hydrodynamic modes due the nonlinear 
term in Eq. (2^). We treat that term as a perturbation, 
i.e. we formally multiply it by a coupling constant 7 
(whose physical value is unity), and calculate its effect to 
flrst order in 7. 

We flrst deflne the velocity autocorrelation function. 



Co.p{k,t) = (u„(k,i) u;3(-k,0)) 



(2.5) 



with k the spatial Fourier transform variable. An equa- 
tion for C can be obtained by Fourier transforming Eq. 
( p.4[ ), multiplying by u^(— k, 0), and averaging over the 
noise while keeping in mind that the noise is uncorrelated 
with the initial fluid velocity. In the case of an incom- 
pressible fluid we need to consider only the transverse- 
velocity correlation function, C±. This is easily done 
by multiplying with unit vectors, fc^'' (i — 1, . . . , d — 1), 
that are perpendicular to k, which eliminates the pres- 
sure term. We obtain 

{dt + jyk^) Ci(k, t) = -nfc^fcW fcf^'Jj ^(u^(k - q, t) 



xuaiq, t) u/3(-k,0)) 



(2.6) 



wher e we have used the incompressibility condition, Eq. 
( ^.31) , to write all gradients as external ones. To zemLh. 
order in 7 we flnd, with the help of the f-sum ruleEj'L2l 
Ci(k,t = 0) ^fceT/p, 



Ci(k,t) = (fcBT/p)e-'"'* + 0(7) 



(2.7) 



This is the standard result obtained from linearized hy- 
drodynamics, which predicts exponential decay for k 7^ 0. 
To calculate corrections due to the nonlinearity, we need 
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an e qua tion for the three-point correlation function in 
Eq. (2_^). The simplest way to do this is to use the time 
translational invariance property of this correlation func- 
tion to put th e tim e depeadence in the last velocity, and 
then use Eq. (2.4) againE^I The result is an equation for 
the three-point function in terms of a four-point one that 
is analogous to Eq. (2.7). To solve this equation, we note 
that, due to the velocity being odd under time reversal, 
the equal-time three-point correlation function vanishes. 
By means of a Laplace transform, one can therefore ex- 
press the three-point function as a product (in freq uency 
space) of the zeroth order result for C_l, Eq. (p?7|), and 
the four-point function. To leading (i.e., zeroth) order 
in 7 the latter factorizes into products of velocity auto- 
correlation functions. Upon transforming back into time 
space, and to quadratic order in the coupling constant 7, 
we thus obtainEJ 

(dt + iyk^) C^(k,t)+ f dr I](k,i-r)C_L(k,r) = , 
"'0 

(2.8a) 

with 



X C^j, (k - p, i) + Cau (p, t) (k - p, i) 

+0(7') 



(2.8b) 



The self-energy E is proportional to k^, and thus pro- 
vides a renormalization of the bare viscosity v. For later 
reference we mention that diagrammatically the above 
contribution to E can be represente d by the one-loop 
diagram shown in Fig. |^. Equation ( 2.8bD , and its rep- 
resentation by Fig. ^ illustrate the meaning of the term 
'mode-coupling' as explained in the Introduction. 




FIG. 2. Leading contribution to the self energy E. The 
internal propagators are given by the transverse velocity au- 
tocorrelation function. 

The source of the LTTs is now evident. In our model 
of an incompressible fluid, only the transverse component 
of Cap is nonzero. 



Ca/3(p,t) = (5q/3 -PaPp) C±{p,t) 



(2.9) 



Putting 7 = 1, and defining Si'{t) = limk-,o E(k, i)/k'^, 
one obtains for asymptotically long times 



p \d(d + 2)J [^Trvty/ 



1 



(2.10) 



This is the well-known contribution of the_.transverse- 
velocity modes to the LLT of the viscosity.L3 In a com- 
pressible fluid, a similar process coupling two longitudi- 
nal modes also contributes to the leading LTT. The other 
transport coefficients, e.g., A and ^, also show LTTs pro- 
portional to 1/t'*/^, and all of them have alscJess leading 
LTTs proportional to l/t(''+i)/2 or weaker Jjla 

The correction to the kinematic viscosity is given by 
the time integral over 5v{t)^ and for the frequency- 
dependent kinematic viscosity the Fourier transform of 
5v{t) implies a nonanalyticity at zero frequency. More 
generally, for the frequency or wavenumber-dependent 
transport coefficients, the LTTs imply the asymptotic 
forms 



M^)/m(o)-i-c;^c.('^-2)/2 + 

Mk)M0) = l-6^|k 



d-2 



(2.11a) 
(2.11b) 



where /i represents any of the transport coefficients, ^, 
A, or Z), in a fluid. The prefactors and are positive. 

Of particular interest for the critical dynamics near 
the liquid-gas critical point is a contribution to 5X from 
the coupling of the transverse- velocity fluctuations to the 
entropy fluctuations. Before performing the wavenumber 
integralpthis particular contribution, which we denote by 

a^.(k,t) = ^^;,(p)^(k.pW)2e-(^P-+^-P^.)* . 

P P i 

(2.12) 

Here Dt = \/ pcp is the thermal diffusivity in terms of 
Cp, the specific heat at constant pressure, x is the order 
parameter susceptibility for the phase transition, where 
we have anticipated that near the critical point we will 
need the momentum-dependent x, and p± = p ± k/2. 
Setting k = and carrying out the momentum integral 
leads to a t~'^/'^ LTT. The correction to the thermal con- 
ductivity is obtained by integrating SXit) over all times, 
as in the case of the kinematic viscosity. 

These results imply that for d < 2, conventional hydro- 
dynamics does not exist. Indeed, it is now known that for 
these dimensions the hydrodynamic equations are nonlo- 
cal in space and time. For. a discussion of this topic we 
refer the reader elsew here . B 

By examining Eq. ( 2.81 ) or ( 2.12| ) one easily identifies 
a mechanism by which the LTT effects can become even 
stronger. Consider a system with long-range static cor- 
relations, for instance due to Goldstone modes, or due to 
the vicinity of a continuous phase transition. In either 



case, some susceptibilities, e.g. the x i^i Eq. ( 2.12 ), be- 
come long-ranged, amplifying the LTT effect. In the next 
section we discuss the realization of this scenario in the 
vicinity of a phase transition. We note that in certain liq- 
uid crystal phases some susceptibilities behave as 1/k^, 
amplifying the LTT effect and causing a braskkdown of 
local hydrodynamics for all dimensions d < 4.EJ 
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We make one last point concerning the LTTs. So far 
we have stressed the leading tails that decay like t~'^/^, 
but their are numerous subleading LTTs as well. Most of 
them are uninteresting, but one becomes important near 
the critical point, via the mecha nism discussed in the last 
paragraph. According to Eq. ( 2.12| ), a central quantity 
for determining the critical contribution to A is the shear 
viscosity 77 (which enters v). It turns out that the con- 
tribution to 77 that is dominant near the critical point 
is a subleading LTT away from criticality. It involveSr-a, 
coupling of two heat or entropy modes and is given byO 



1 



X(P) X(k - 



(2.13) 



with A a constant. Away from the critical point, this 
LTT decays as so in fact it is a next-to-next 

leading LTT. Nevertheless, it is the dominant mode- 
coupling contribution to 77 near the critical poi nt be cause 
of the two factors of x in the numerator of Eq. ( 2.13 ). We 
discuss this enhancement of LTTs by critical susceptibil- 
ities further in the next subsection. 



C. Critical dynamics 

Near continuous phase transitions, fluctuations grow 
and ultimately diverge at the critical point. For the 
liquid-gas critical point the order parameter is the dif- 
ference between the density and the critical density, and 
the divergent fluctuations are the density fluctuations as 
described by the d ensity susceptibility. The susceptibil- 
ity X in Eq. ( 2.12| ) is proportional to this divergent sus- 
ceptibility, which in the Ornstein-Zernike approximation 
behaves as 



X(P) oc 



1 



(2.14) 



with p measured in suitable units, and r the dimension- 
less distance from the critical point. Carrying out the 
time integral, the leading singular contribution to t 
static, wavenumber-dependent thermal conductivity isl 



6X(k) 



(2.15) 



Using Eq. ( |2.14| ) in this equation we see that the homoge- 
neous thermal conductivity is infinite at the critical point 
for all d < 4, diverging as |r|~('*~'''/^. This is result of 
the amplification of the LTT by the critical fl uctu ations. 

By the same mechanism, we see that Eq. ( 2.14 ) leads 
to a logarithmically sin gular contribution to Srj, if we take 
into account that Eq. ( ^.15 ) implies that at the critical 
point, Dxik) ~ |k|'*-2. 



The above are one-loop calculations that use the 
Ornstein-Zernike susceptibility. To go beyond these ap- 
proximations one needs to (1) use the correct scaling 
form for the susceptibility, and[-(2) use either a self- 
consistent mode-mi.ipling theory,Ej or a renormalization 
group approachollj [to improve on the one- loop approx- 
imation. The resultEil is that the thermal conductivity 
diverges like 



(2.16a) 



Acx|r|-0-57 , 

and that the thermal diffusivity vanishes like 

Dt cx |r|"-s^ . (2.16b) 

This is in good agreement with experimental results, as 
shown in Fig. ^ which compares experimental and theo- 
retical results for the thermal diffusivity of carbon dioxide 
in the critical region. 
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FIG. 3. The thermal diffusivity Dt = ^/pcp of carbon 
dioxide in the critical region as a function of density at various 
temeratures {Tc — 304.12 K). The symbols indicate experi- 
mental data for Dt measured directly, and for X/pCp deduced 
from thermal-conductivity data. The solid curves represent 
values calculated from the mode-coupling theory. From Ref. 



D. Nonequilibrium effects 

We now turn to a fluid in a nonequilibrium steady 
state. It has been known for some time that these sys- 
tems in general exhibit GSI in spatial correlations, and 
that the order parameters for phase transitions in such 
systems couple to these long-range correlations. The first 
study of a phase transition of this type considered phase 
separation in a binary liquid under shearO The spatial 
correlations responsible for the GSI are closely related 
to those that lead to the LTTs in equilibrium time cor- 
relation functions. We will consider a fluid in a steady, 
spatially uniform temperature gradient VT, but far away 
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from any convective instability. Further, we use a num- 
ber of approximations that enable us to focus on the most 
interesting effects of such a temperature gradient. For a 
justification of this procedure, as well as the underlying 
details, we refer to Rcf. 

Focusing on the coupling between fluctuations of the 
transverse fluid velocity uj_ and temperature fluctuations 
ST, and neglecting the nonlinearity in Eq. (2.11;), Eqs. 
(2.1) can be written 



dt + (u • V)T 



1 



I^V^M±,„ + — (V/3Pa/3)_L , (2.17a) 



Po 



DtS7^ST 



PoToCp 



(V-q) , (2.17b) 



where po and Tq indicate average values. These bilinear 
equations can be solved by means of Fourier and Laplace 
transformations. Focusing on sta,tic or equal-time corre- 
lations, one finds, for example,La 



(|<5p(k)p)=pfcBT 



dp 

dp 



pkBT{aTkj_ ■ \JTf 



((k^ •g(k))(5p(-k)) -pfceTaT 



kx ■ vr 



with, 



aT 



-1 
P 



dp 
dT 



(2.18a) 
(2.18b) 

(2.18c) 



the thermal expansion coefficient. 

There a re seve ral remarkable aspects of these results. 
First, Eq. (2.18a) for the density correlations implies that 
the first term, which also exists in equilibrium, is delta- 
function correlated in real space, while the second term 
decays like const. — Ixlicuthree-dimensions, with |x| the 
distance in real spaccEjO Equation ( 2.18b| ) shows that 
the transverse momentum-density correlation functions 
decays as l/|x| in three-dimensions. Both of these re- 
sults show that spatial correlations in a NESS exhibit 
GSL Second, the right-hand side of Eq. (|l8^) IS essen- 
tially the integrand in Eq. (2.15), which shows the close 
connection between the LTTs, singularities in transport 
coefficients near continuous phase transitions, and the 
GSI of spatial correlations in a NESS. 

The nonequilibrium fluctuations arise from a cou- 
pling between the temperature fluctuations perpendic- 
ular to the temperature gradient and the transverse- 
momentum (vis c o u s ) fluctuations parallel to the temper- 
ature gradient .eSEjo The amplitudes of these fluctua- 
tions can be measured by small-angle light-scattering ex- 
periments. Specifically, the dynamical structure factor 
S'pp(k, i) = (5p(k, t) p(— k, 0)), which is proportional to 
the scattering cross section, in a NESS has the formH 



Spp(k,t) = So (l + ^T)e 



A, 



(2.19a) 



where Sq is the structure factor in equilibrium, and 



Cpv/Dr (k± • VT)2 



Av = 



(ki • VT)2 



T(i/2 - £)2 ) 



k4 



(2.19b) 
(2.19c) 



For t — Q one rec overs t he equal-time density correla- 
tion function, Eq. ( 2.18a ). The amplitudes Ax and A^ 
are proportional to (Vr)f/k'*, which has been verified 
by experiments, see Fig. ]^ The agreement, with no ad- 
justable parameters, is excellent. Notice that the am- 
plitude of the temperature fluctuations is enhanced by 
a factor of a hundred compared to the scattering by an 
equilibrium fluid. 




k 



FIG. 4. Amplitudes At and of the nonequilibrium tem- 
perature and transverse-momentum (viscous) fluctuations in 
liquid hexane at 25° C as a function of (Vr)Vk''. The sym- 
bols indicate experimenta l dat a. The solid lines represent the 
values predicted by Eqs. (2.19). From Ref. 



In a liquid mixture or a polymer solution, a tem- 
perature gradient induces a concentration gradient due 
to the Soret effect. Nonequilibrium concentration fluc- 
tuations then arise from a coupling between the con- 
centration fluctuations and the transverse-momentum 
fluctuations.Ha Just as in the case of the nonequilib- 
rium temperature and viscous fluctuations, the ampli- 
tude of the nonequilibrium fluctuations will also be pro- 
portional to (VT)^/k^. The presence of such long-range 
nonequilibrium concentration fluctuations has also been 
confirmed experimentally, see Fig. ^. 

One final point concerns the question what these 
generic long-range spatial correlations would do to criti- 
cal behavior. This is a natural question, as usually long- 
range interactions, or long-range correlations, fundamen- 
tally modify important aspects of any phase transition. 
This questions has been addressed elsewhere in the con- 
text of phase transition in nonequilibrium, or driven. 
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Vf I k' (xW^" cm) 

FIG. 5. Amplitude Ac of nonequilibrium concentration 
fluctuations in polystyrene and toluene solutions at 25° C as 
a function of (Vr)^/k*. The symbols indicate experimental 
data at polymer mass fractions w = 4.00% (solid symbols) 
and w — 0.50% (open symbols), respectively. Different sym- 
bol shapes correspond to different values of the wave number 
k. The solid lines are linear fits to the exjwimental data in 
agreement with the theoretical predictions Mj 

classical systenis.@ In the remainder of this paper we 
study the analogous question in the context of quantum, 
or zero-temperature, phase transitions in equilibrium sys- 
tems. 



III. GENERIC SCALE INVARIANCE AND 
PHASE TRANSITIONS IN QUANTUM SYSTEMS 

A. Field theory 

We now turn to a quantum fluid, specifically, to sys- 
tems of interacting electrons in solids. Since we will 
again be interested in effects at long wavelengths and 
long times, the ionic lattice and effects resulting from it 
will be irrelevant for our purposes, and we therefore adopt 
a simple 'jellium' model of free electrons interacting via 
the Coulomb interaction in the presence of a charge neu- 
tralizing homogeneous background. There are various 
ways to theoretically deal with this many-body problem, 
but the most powerful one for our purposes is a field- 
theoretic forpuilation of quantum statistical mechanics. 
This methodllj starts by writing the partition function in 
terms of a functional integral, 



Z = 



V>(0) = -i/'(l/fcBT) 



(3.1) 



Here D is a functional integration measure with respect 
to fermionic, i.e., Grassmann valued, fields ip and ^ that 
are defined on a real-space/imaginary-time manifold with 
the time sector spanning the interval [Ojl/ksT]. The 
action S for our model system reads 

dx ^i^a{x) [--a^ + V V2m + ^] V'aW + S'int . 

(3.2) 



Here x = (x, r) with x denoting position, r imaginary 

time, and a is the spin index. J dx = J dx j^^^^"^ dr, 
m is the electron mass, and fi is the chemical potential. 
At T = 0, /i = kp/2m with k-p the Fermi wavenumber. 
Here and in what follows we use units such that h = 1. 
At the most basic level of the theory, S'int represents the 
Coulomb potential, but it often is advantageous to inte- 
grate out certain degrees of-freedom to arrive at effective, 
short-ranged interactions.EZI For our purposes we do not 
need to specify the precise form of S'int, it will suffice to 
postulate that the ground state of the interacting system 
is a Fermi liquid. 

It is useful to go to a Fourier representation with 
wavevectors k and fermionic Matsubara frequencies w„ = 
27rT(n-|- 1/2), 



z(kx 



(3.3a) 
(3.3b) 



From Eq. (3.2) it follows that, in a free-electron system. 



single-particle excitations about the Fermi surface are 
soft with a linear dispersion. 



z 



ILUn -Vpq 



(3.4) 



with Z = I, q = — kp, and v-p ~ fcp/m the Fermi 
velocity. Here (...)= / Dlip, ^j] . . . e'^/Z comprises both 
a quantum mechanical and a therm ody namic average. 
Fermi liquid theory shows that Eq. ( |3.4| ) remains valid 
in the presence of S'int, only the values ot-|Z and vp 
are changed compared to a free electron gas.H Similarly, 
there are soft two-particle excitations. For instance, the 
susceptibility 

k,p n' ,m' 

>^ {-ipn' .a' {P - q)i^m' {P))) (3.5) 

diverges for q = 0, nm < 0, and i^n-m = ^ ^ 
like l/|r2„_„i|. At nonzero external wavenumbers |q|, the 
wavenumber scales like the frequency, so the dispersion 
of the soft modes is again linear. 

In many solid-state electronic systems, quenched disor- 
der is present in the form of impurities or lattice defects. 
This has important consequences for the transport and 
thermodynamic properties of the system, and we model 
it by means of a term in the action 



S'dis = / dx tt(x) y ^-0g.(a:) V'cr(x) 



(3.6a) 



Here u(x) is a static, random potential governed by a dis- 
tribution P[u]. For simplicity, we take u(x) to be delta- 
correlated and Gaussian distributed with the second mo- 
ment of P given by 
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{M(x)M(y)}dis = 



1 



TrA^FTrol 



(5(x- 



(3.6b) 



Here {. . .}dis denotes the disorder average, r^d is the elas- 
tic relaxation time, and A^f is the density of states at the 
Fermi surface. I n the presence of S'dis7 the single-particle 
excitations, Eq. ( |3.4| ), are obviously massive with a mass 
proportional to l/rrd. However, two-particle excitations 
are still soft. In contrast to the clean case, their dis- 
persion relation is diffusive. For instance, the disorder 
averaged analog of the-two-particle propagator given in 
Eq. (|3.5D has the fornO 



E E {((Vw(k + q)V^„.,.(k)) 



k,p n',m' 
X 



(V'n',<T'(P - q)V^m',<T'(p)))| OC — 5- 

J dis |ii,i_„i| - 



Z3q2 ' 
(3.7) 



where the diffusion coefficient D is proportional to Tioi. In 
both the clean and the disordered cases, the two-particle 
excitations have the same form for different combinations 
of the spin indices. 

Notice that the softness of the two-particle excitations 
is not related to any conservation law, i.e., the combi- 
nations of fermion fields do not correspond to a density 
or, in the clean case, a current. Rather, these excitations 
are the Goldstone modes resulting from the spontaneous 
breaking of the symmetry between positive and negative 
Matsubara frequencies. The order parameter that be- 
longs to this symmetry is the quantity 



'= lim (V'„.cr(x) o.(x)) 

- lim (?/5„.c,(x)?/'n.<T(x)) 



(3.8) 



which is the single-particle spectral function, or the dif- 
ference between the retarded and advanced Green func- 
tions. Since the causal Green function has a cut on the 
real axis, Q is nonzero as long as the density of states at 
the Fermi surface is nonzero. This connection between 
the single-particle spectral function and the soft mode 
specteum was discovered by Wegner for the disordered 
casCjEa and has been elaborated on in Refs. |50| , |5l| , ^ . 



B. Long-time tails and spatial GSI in equilibrium 

1. Long-time tails, a.k.a. weak-localization effects 

At zero temperature, a clean electron system has an 
infinite diffusivity or conductivity, so we start our dis- 
cussion with the disordered case. One expects the mode- 
coupling argument given in Sec. || to still apply, with 
the transverse-velocity mo des of Eqs. replaced by 

the diffusive modes of Eq. (3.7). Since the dispersion re- 
lation is the same in both the classical fluid case and the 
disordered electron case, we thus expect a LTT in the real 



part of the c onduc tivity that takes the same functional 
form as Eq. ( ^.lla ), viz. 



(3.9) 



This expectation is indeed borne^out by an explicit cal- 
culation for both noninteracting^ and interactingo elec- 
trons. It turns out, however, that the coefficient cJJ is 
negative, in contrast to the classical case. This difference 
in the sign of the LTT is due to the fact that the scat- 
terers that lead to a flnite conductivity in the quantum 
case are static, while in the classical fluid they are the 
moving fluid particles themselves. The strength of the 
LTT, on the other hand, is the same in both cases due to 
the strong similarity of the respective soft mode spectra. 

In this context it is interesting to point out that an- 
other, and in some sense closer, classical analog of the 
quenched disordered electron fluid is the classical Lorentz 
gas, which consists of noninteraeting classical particles 
moving between static scatterersBil The diffusivity in this 
case also shows a LTT with the same sign as in the 
quenched disordered quantum case (and for the same rea- 



son), but the strength of the LTT is weaker, viz. w 



d/2 

in frequency space or time space. The rea- 

son is that in the classical Lorentz gas there is no analog 
of the spontaneously broken symmetry mentioned above, 
and the only soft mode (wliich is diffusive) is due to par- 
ticle number conservation.!^ There are thus many fewer 
soft modes than in either a classical fluid or a disordered 
electron system, which leads to a weaker mode-coupling 
effect and, hence, to a weaker LTT. Indeed, at nonzero 
temperature the uj^'^~'^y^ dependence in Eq. (3.9) gets 
transformed into a T^'*^^^/^ behavior, and in the limit 
(jj ^ at fixed T > one recovers the result for the clas- 
sical Lorentz gas. This example illustrates that for the 
purpose of comparing LTT effects in different systems 
the crucial criterion is how similar the soft-mode spec- 
tra are. In this respect the disordered electron system, 
even without electron-electron interactions, is closer to 
the classical ffuid than to its classical limit, the classical 
Lorentz gas. 

The fact that a nonzero temperature cuts off the LTT 
singularity provides the most convenient way to experi- 
mentally observe the phenomenon in electronic systems, 
since the dynamical conductivity is hard to measure. Fig. 
^ shows the VT dependence of = 0, T) in bulk Si:B. 
In 2-d systems, r ealiz ed by thin metallic films, the zero 
exponent in Eq. ( |3.9| ) translates into a logarithmic fre- 
quency or temperature dependence within a large dy- 
namical range, like in the classical case, Eq. ( p^ . An 
experimental example of this phenomenon is shown in 
Fig. 0. 

These LTTs in disordered electron systems are known 
as 'weak- localization effects'. This is because they can 
be considered precursors of the 'strong localization' that 
occurs in systems with very strong disorder, which un- 
dergo a phase transition froin_a metal to an insulator 
as a function of the disorder.^ In two-dimensions the 
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FIG. 6. Conductivity data from Ref. ^ showing the static 
conductivity of ten Si:B samples in a magnetic field plotted 
against y/T at low temperatures. From Ref. p3. 
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FIG. 7. Resistance data from Ref. for a thin PdAu film. 
The resistance R, normalized to Ro ~ R{T = IK), is plotted 
versus logT. From Ref. B8. 



weak-localization effects are strong enough_.to prevent 
the formation of a true metal altogether,E3 and as a 
function of increasing disorder the system undergoes a 
crossover from a weak insulator to a strong one. Systems 
in d < 2 are always localized, or insulating, irrespective 
of the strength of the disorder, at least for noninteracting 
electrons. 



2. Generic scale invanance vn equilibnum 

In quantum statistical mechanics, as opposed to the 
classical theory, there is an intrinsic coupling between 
statics and dynamics. This is manifested by the basic 
fermionic fields ip and "0 being functions of both space 
and (imaginary) time. Consequently, one expects effects 
similar to the LTTs in time correlation functions to ap- 



pear in static correlation functions, even in equilibrium. 
As an example, let us consider the wavenumber depen- 
dent static spin susceptibility Xs in a disordered electron 
system. Given that the wa venu mber scales like the square 
root of the frequency, Eq. (3.7), it is natural to guess that 
Xs at small wavenumbers has the form 



Xs(q) = Co - cd_2|q 



d-2 



O(q') 



(3.10) 



The coefficients cq and Cd~2 are both expected to be 
positive. This is because the quenched disorder slows 
the electrons down, which leads to an increased effective 
electron-electron interaction. Consequently, the disorder 
enhances cq compared to the value of the homogeneous 
spin susceptibility in a clean system, and the suscepti- 
bility decreases with increasing wavenumber. This is the 
sam e effect that leads to a positive coefficient in Elq, 
( |3.9D . This result is confirmed by explicit calculations .e2I 
This nonanalytic wavenumber dependence, which for di- 
mensions 2 < c? < 4 is the dominant one, corresponds 
to a l/r^'^''"^) decay in real space. We see that the cou- 
pling between statics and dynamics produce GSI in quan- 
tum systems even in equilibrium. Note that this effect is 
much weaker t han t he GSI seen in non-eq ui libriu m clas- 
sical (see Sec. IIP ) or quantum (see Sec. HID below) 
systems. 

Detailed calculations show that the origin of the GSI 
can be traced, as the above argument suggests, to the 
same Goldstone modes that lead to the LTTs and were 
discusssed in the previous subsection. The relevant con- 
tributiQiL|to Xs can be schematically represented by the 
integrate 



dpp' 



|q| 



duj 



(^ + T+p2)3 



(3.11) 



with A a microscopic wavenumber. Equation (3.11) 
demonstrates the coupling between statics and dynamics 
mentioned above. Notice that for power-counting pur- 
poses the integrand is again a product of t wo dif fusive 
modes. At T = 0, the integral yields Eq. ( |3.10 ), and 
at T > the temperat ure cut s off the leading singular- 
ity as discussed in Sec. [II B 1 above. In contrast to the 



LTT of the previous subsection, the above results hold 
for interacting electron systems only. For noninteracting 
electrons, the absence of frequency mixing prevents the 
coupling of the static spin density fluctuations to the soft 
modes, and Xs is analytic at zero wavenumber. 

One might expect all static correlation functions to 
display such nonanalytic wavenumber dependences, but 
this turns out not to be true. For instance, the parti- 
cle number density susceptibility does not have a leading 
singularity analogous to that in Xs , and neither does any 
other spin-singlet particle- hole susceptibility like, e.g., 
the number density current susceptibility. This is be- 
cause these observables couple less strongly to the soft 
modes than the spin density. In Ref. |6l| general crite- 
ria have been developed that allow to determine which 
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susceptibilities show GSI due to soft modes, and which 
do not. This reference also discusses limitations of the 
equivalency between statics and dynamics in quantum 
statistical mechanics. For instance, it is remarkable that 
the conductivity in a disord ere d electron systems shows 
an uj^''-^'^)/'^ LTT, see Eq. (^), but the corresponding 
static current susceptibility does not show an analogous 
|q|'^~2 wavenumber dependence. This is because a finite 
frequency breaks the symm etry in frequency space that 



we mentioned in Sec. Ill A. As a result, dynamical cur- 



rent fluctuations couple more strongly to the diffusive 
soft modes than static ones. 

We now turn to the clean case. As we have seen, the 
only difference compared to the disordered case is the 
dispersion relation of the soft modes, which now is bal- 
listic rather than diffusive, i.e., wavenumbers scale like 
frequencies, which leads to weaker LTT and GSI effects. 
Let us aga in consider the static spin susceptibility. From 
Eq. (3.11) we expect that, at zero temperature, there is 
a contribution to Xs that is of the form 



A 
'|q| 

wich leads to 
Xs(q) = 



dpp 



d-l 



duj 



Cd-i|q| 



id-l 



0(q2 



(3.12a) 



(3.12b) 



with positive coefficients Cq .and Cd-i- Explicit calcula- 
tions confirm this result.HS As in the case of the LTTs, 
the signs of the nonanalyticities in the clean and disor- 
dered cases, respectively, are different. A simple physical 
explanation in the present case is that the nonanalyt- 
icity is produced by fluctuation effects that weaken the 
tendency towards ferromagnetism. Consequently, they 
decrease the value of Xs(q = 0)i and with increasing 
wavenumber the susceptibility increases. 



C. Quantum phase transitions 



lie 



In Sec 
ical behavior 



we have seen how LTTs affect the crit- 
at a classical phase transition. Here we 
discuss the., analogous problem for a quantum phase 
transition.E3L3 An example is thCpiprromagHctic transi^ 
tion that is observed in, e.g., MnSi,E3 UGe2,E3 or ZrZn2E£l 
at T = as a function of hydrostatic pressure. A crucial 
difference between the quantum and classical cases is that 
in the latter, due to the lack of couphng between statics 
and dynamics, the LTTs affect only the critical dynamics 
and the critical behavior of quantities given by time cor- 
relation functio ns, l ike the thermal conductivity in the 
example in Sec. II C . In the quantum case, however, the 
LTT/ GSI phenomena also influence the critical behavior 
of thermodynamic quantities. The easiest way to see this 
is by deriving an order-parameter or Landau-Ginzburg- 
Wilson (LGW) theory for the phase transitionO 



1. Order-parameter field theory 

Let us consider a phase transition with an order- 
parameter field n{x) — n{'ip(x),'ijj{x)) that is bilinear 
in the fermion fields. For example, in the ferromag- 
netic case, n is a vector field, viz. the spin density 
"iT-ii^) = Z]<T,o-' V'<T(a:)(o-i)(T,o-''0<T'(a;), with at [i = x,y,z) 
the Pauli matrices. In general, n is a rank-m tensor field. 
In order for a phase transition to a phase with a nonva- 
nishing expectation value of n to occur, the interaction 
term 5int in the action must contain an interaction be- 
tween the order-parameter modes. This term, which we 
denote by S^f , reads schematically 



5Sr=r / dxn'ix) 



(3.13) 



with r a coupling constant. For simplicity, here and in 
what follows wc suppress both the spin and tensor labels, 
and use an obvious symbolic notation. The full action we 
write as 



5* — Sq + S; 



OP 



(3.14) 



with the action S'o, which we will refer to as the 'reference 
ensemble', containing alL contributions other than . 

We now follow HertaS3 in deriving an order-parameter 
field theory for the phase transition under considera- 
tion. To this end, we first decouple by means of 
a Hubbard- Stratonovich transformationEj Denoting the 
Hubbard- Stratonovich field by M, we write the partition 
function 

Z = f D[i>,^] e^I'^^'^l 

const.x ID[M] (^^-2r J d.M(.) ni.)^ 

(3.15) 



const. X / D[M] e-*[^^l , 



where (. . .)o denotes an average with the reference en- 
semble action Sq, and $[Af] is the LGW functional. The 
latter reads explicitly 



<^>[M] =T I dx Wix) - In^e-'^/''" M(x)„(.)^^ 



(3.16) 



and can be expanded in powers of M, 

1 



$[M] = dxidx2 M{xi) 



r 



S{xi - X2) 



X^'^\xi - X2) 



M{X2) 



+ Y\ j dxidx2dxz x^^\xi,X2,X'i) 

xM{xi)M{x2)M{x3) + 0{M*) , (3.17a) 
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where we have scaled M with l/^ /2T. The coefficients 



X^'^ in the Landau expansion, Eq. (|3.17a ), are connected 



l-point correlation functions of n{x) in the reference en- 
semble, 



X^^Kxi, ...,xi)^ {n{xi) ■ ■■n{xi))[ 



(3.17b) 



2. Local versus nonlocal LGW theories 

A crucial question now arises concerning the behav- 
ior of the correlation functions defined in Eq. ( 3.17b ). 
Suppose their Fourier transforms are finite in the limit 
of small wavenumbers and frequencies. In lowest order 
in a gradient expansion one can then simply localize the 
cubic and higher-order coefficients. Defining a Fourier 
transform of the Hubbard-Stratonovich field by 



M{q) 



Vt/vJ 



dx e 



M{x) 



(3.18) 



we obtain an LGW functional 

m] = lT. ^(«) [i/r - x^'^ (q) M{~q) 



+ 1t5 II Siqi+q2 + q3 + qi)M{qi)M{q2) 



4! V 



xM{q3)M{q4) 



(3.19) 



Here the four- vector q ~ (q, cij„) comprises wavevector 
and Matsubara frequency, and we have assumed that x'^' 
vanishes at zero frequency and zero wavenumber. Sup- 
pose x^^"^ is an analytic function of the wavenumber. 



X^'Hq) = co + C2Cl' + cn\nn\/\ci\ 



(3.20) 



where the value of m depends on the specific system. 
Then we obtain an ordinary Landau theory with the 
Gaussian part of the LGW functional given by 

'^^^^M] = i ^ M{q) [r + a2 q2 -f ao|f7„|/|qr] 

9 

xAf(-g) . (3.21) 

Here r ~ l/L — x'^^H'l — 0) the bare dista nce from 
the critical point at T = 0. Equation ( ^.21 ) leads to 
mean-field values for the critical exponents v and 7, and 
a dynamical critical exponent z = m + 2. Power count- 
ing shows that the non-Gaussian terms are irrelevant in 
the renormalization-group sense of the word and, hence, 
the exact critical behavior is also mean fieldjiy^e. This 
is the conclusion that was reached by Hertz J2a namely, 
that generically quantum phase transitions have mean- 
field critical behavior. 



As we have seen in Sec. [HB2, however, this scenario 



does not necessarily apply. Equation (3.1C) provides an 
example of a susceptibility that does not have the form 



of Eq. (3.20) due to GSI effects. Moreover, if the nonan- 
alytic wavenumber dependence of the susceptibility x'^^ 
is caused by soft modes that are made massive by an ex- 
ternal field H conjugate to the order parameter (or by a 
nonvanishing average of the order parameter), then the 
singular part of the static, field-dependent susceptibility 
will have the form 



0,H) 



H)y 



(3.22) 



where y > and a; > are exponents that determine the 
nature of the nonanalyticity and the scaling of H with 
the wavenumber, respectively. This is true, for instance, 
in the case of the disordered itinerant ferromagnet, where 
X = 2 and y = {d— 2)/2. Since the higher susceptibilities 
X^^\ X^^h etc., can be obtained from x^^-* by differenti- 
ating with respect to H, this implies that they diverge 
in the limit of zero wavenumbers and frequencies. That 
is, the gradient expansion that leads to the usual LGW 
theory, and hence LGW theory itself, do not exist. In- 
stead, one obtains a field theory where the coefficients 
of the various powers of the field M are singular func- 
tions of wavenumber arui-frequency. In other words, the 
field theory is not local.t2l The above arguments make it 
clear that this is a direct consequence of the GSI in the 
equilibrium quantum system. 

Such nonlocal field theories are difficult to handle and 
not suitable for explicit calculations. For the purpose of 
determining the critical behavior one can try to circum- 
vent this problem by checking whether the field theory 
still allows for a Gaussian critical fixed point with respect 
to which the higher-order terms are irrelevant by power 
counting. This was the approach taken in Ref. 60 for the 



disordered ferromagnetic transition, and it relied strictly 
on power counting by taking the divergence of, e.g., the 
coefficient U4 in Eq. (3.19) into account by assigning it 
a suitable scale dimension. It turns out that this proce- 
dure yields the correct power laws, but misses logarithmic 
corrections to scaling that power counting is not sensitive 
to. In the following subsection we use this example to il- 
lustrate the interplay between GSI and quantum critical 
behavior. 



3. Example 1: The disordered quantum ferromagnetic 
transition 

As an example of the interplay between LTT/GSI ef- 
fects and quantum critical phenomena, we now consider 
the quantum critical behavior at the ferromagnetic tran- 
sition in quenched disordered itinerant electron systems 
in d-dimensions. As we mentioned in the last subsection, 
in that case the two-point susceptibility x^^^ has the form 



X^^^g) = Co + Crf_2|q|'' ^ + C2q^ + cn|f^„|/q' 



(3.23) 



This leads to a Gaussian part of the LGW functional 
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+an|f^„|/|q|' M{-q) . (3.24) 



Furthermore, the fi eld-d ependent susceptibihty has the 
form given by Eq. ( |3.22| ) wit h x = 2 and y = (d - 2)/2. 
This imphes that 7/4 in Eq. ( |3.19D , and aU of the coeffi- 
cients of higher-order terms in the Landau expansion, do 
not exist in the hmit of small wavenumbers and frequen- 
cies. For instance, U4 behaves in this limit like 



U4(|q| ->0,1]„ = 0) 



(d-6) , 



|d-6 



,(0) 



(3.25) 



with u^^ and u]^' finite numbers. 

Let us now look for a Gaussian fixed poin t that de- 
scribes the critical behavior. From Eq. ( 3.24 ) one reads 
off the Gaussian values of the correlation length expo- 
nent v, the susceptibility exponents 7 and 77, and the 
dynamical exponent z as 

i/=l/{d-2) , 7=1 , rj^A-d , z^d , 

(3.26a) 



,(0) 



for 2 < d < 4, and 



1/2 



7 = 1 



?7 = 



z = 4 



(3.26b) 



for d > 4. In d = 4 one finds mean-field exponents with 
logarithmic corrections to scaling. Power counting sug- 
gests that the non-Gaussian terms in the action are ir- 
relevant with respect to this fixed point, so that these 
exponents represent the exact critical behavior. How- 
ever, for the remaining exponents, the order-parameter 
exponents (3 and 5, one needs to take into account these 
terms, as is the case in ordinary Landau theory. From Eq. 
( 3.22 ) with a; = 2 we see that the average order param- 
eter m, which scales like H , scales like the wavenumber 
squared. The coefficient U4, Eq. ( 3.2£ ), thus scales like 



-d/2 



which suggests an equation of state 



= H 



with u, w > 0. From this we obtain 

^ = 2/(d-2) , 5^ d/2 
for 2 < d < 6, and 

/3 = l/2 , 5 = i , 



(3.27) 



(3.28a) 



(3.28b) 



for d > 6. In d = 6 there are again logarithmic correc- 
tions to scaling. A more sophisticated analysisE^I confirms 
these results, which illustrate the strong influence of the 
GSI on the c ritical behavior: The exponents given in Eqs. 
( 3.26a ) and ( 3.28a ) for the experimentally most interest- 
ing case d = 3 are very different from the mean-fijeki 
exponents one obtains if one neglects the GSI effects.E3 



While the above simple power-counting arguments pro- 
duce the correct exponents of the exact critical behav- 
ior, it turns out that the simple Gaussian fixed point 
presented above is marginally unstable. This has been 
traced to the existence of two time scales, viz. the criti- 
cal one with dynamical exponent 2;, and the diffusive one 
with dynamical exponent 2. The existence of the latter 
has been obscured in the derivation of the LGW theory, 
which integrates out the diffusive modes. A more care- 
ful analysis, which keeps all of the soft modes explicitly 
and on equal footing, shows that the actual critical fixed 
point is not Gaussian, but has the same exponents as 
given above, with multiplicative logarithmic corrections 
to scaling.C2l A convenient way to account for these cor- 
rections is to write the critical behavior as power laws, 
with scale-dependent critical exponents. For instance, 
the correlation length ^ as a function of the distance r 
from criticality has a log-log normal correction to the 
simple power law. 



(3.29a) 



where g{x) is a function whose leading behavior at large 
arguments is 



']{x —t 00) = const. X e 



[ln(c(d)x)]^/2 1n(d/2) 



(3.29b) 



with c(d) a dimensionality dependent constant. This be- 
havior can be represented by writing ^ c>c r^^/^ with a 
scale dependent exponent z/. 



l/i^ = d-2 + ln,g(&)/ln& 



(3.30a) 



with h an arbitrary renormalization group length scale 
factor. For instance, to translate the dependence of h 
into an r-dependence, one needs to substitute b — r^'^ . 
For the other exponents, the analysis of Ref. [70| yields 



z = 2(5 = d + ln.g(6)/ln6 
ry = 4 — d — In g(6)/ In 6 



f3 = 2v 



(3.30b) 
(3.30c) 
(3.30d) 



Equations ( |3.29D and ( |3.30D are valid for 2 < d < 4. The 
exponent 7 has no logarit hic c orrections, and is 7 = 1 
exactly, as stated in Eqs. ( 3.26 ). 

The above detailed predictions have yet to be tested 
exp eriment ally. 



4-. Example 2: The clean quantum ferromagnetic transition 

As a second example we briefly consider the case of 
a clean ferromagnet. In this case, the coefficient U4 of 
the quartic term in the Landau expansion diverges like 
l/|q|''~^, and the average order parameter m scales like 
the wavenumber. Furthermore, as we have mentioned in 



Sec. [II B 2 , the sign of the nonanalytic term is opposite 
to that in the disordered case. This leads to an equation 
of state 
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rm — vm + u m 



for general d, and 



rm + v In + ' 



H 



= H 



(3.31a) 



(3.31b) 

for d = 3, with u,v> O.EI 

These considerations suggest that, generically, the 
quantum ferromagnetic transition of clean itinerant elec- 
trons js-iof first order. Indeed, this-4^ what is observed in 
MnSijEa see Fig. ^, and in UGe2,E3 both systems where 
the ferromagnetic transition can be triggered at very low 
temperatures by applying hydrostatic pressure. Both a 
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where (. . .) denotes a nonequilibrium thermal average. 
The ca lculation shows that the Fourier transform of Eq. 
( |3.32aD , Ci(k), behaves like 



zk • V/i 
2nDk^ 



(3.32b) 



In real space, this corresponds to a decay proportional 
to l/jxl"^"^. This is the same result as the one in the 
corresponding classical Lorentz gas, since the additional 
soft modes in the quantum system do not couple to Ci . 

It is also interesting to consider the electronic structure 
factor, 



C2(x,y) ={((5n(x) My)>}c 



(3.33) 



Let us consider the nonequilibrium part of C2. As one 
would expect, the calculation yields a result that is anal- 
ogous to the one in the classical fluid, Eq. ( ^.18a ), viz. 



C2(k) 



rcl 



6dTT{Dk- 



2^2 



25(V/i)2 - 12(k- V^i)^ 



(3.34) 



As in the classical case, this corresponds to a decay in 
real space like const. — |x| in three-dimensions. For a 
classical Lorentz gas one finds instead 



C2(k) oc (VAi)Vk2 



(3.35) 



The weaker singularity in this case is due to the fact that 
the cla ssical m odel has fewer soft modes, as we discussed 
in Sec. IlIIBl. 



FIG. 8. Phase diagram of MnSi. The insets show the be- 
havior of the susceptibility close to the transition as reported 
in Ref. ^. From Ref. Q 



nonzero temperature and nonzero disorder cut off the sin- 
gularity in Eq. ( 3.31b| ), and disorder, o f cou rse, induces 
the stronger singularity shown in Eq. ( 3.27 ). For small 
values of both temperature and disorder, the competi- 
tion between these terms leads to an interesting structure 
of the phase diagram, with tricritical points and critical 
endpoints. This has been discussed in detail in Ref. 71 



D. Nonequilibrium effects 

Very recently, spatial correlations of density fluctua- 
tions have been studied in noninteractingjiisordered elec- 
tronic system that are not in equilibrium.L3 For the model 
defined by Eqs. (3.2) and (3.6), without the electron- 
electron interaction term Sint but in the presence of a 
chemical-potential gradient V/x, this reference calculated 
an electron density-scatterer density correlation function 
defined by 



Ci(x,y) = {(5n(x))w(y)U 



(3.32a) 



IV. CONCLUSIONS 

In this paper we have emphasized that the soft modes 
which always exist in many-body systems due to either 
conservation laws or broken continuous symmetries via 
Goldstone's theorem, generically couple to the correla- 
tion functions relevant for both scattering and transport 
experiments, and lead to power-law correlations in space 
and time, in the entire phase diagram. That is, correla- 
tion functions usually exhibit generic scale invariance. 

In the classical case we have seen that these GSI ef- 
fects can get amplified either in a phase with Goldstone 
modes, or near continuous phase transitions. In both 
cases the crucial point is that the coefficient of a corre- 
lation function exhibiting GSI is proportional to a static 
susceptibility that diverges either in an entire Goldstone 
phase, or at a special critical point. 

The quantum case is of particular interest because of 
the inherent coupling between statics and dynamics in 
zero-temperature systems. Generically, this implies that 
systems that exhibit GSI in the time domain do so in 
space as well. This has been contrasted to classical sys- 
tems, where the equilibrium fluid provides an example of 
a system that exhibits GSI in time correlations functions. 
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i.e., that has LTTs, but exponentially decaying correla- 
tions in space. This coupling also leads to the important 
conclusion that many quantum phase transitions are in a 
non-mean field universality class governed by long-ranged 
interactions. 

Among the many experimental consequences of the 
particular type of GSI discussed here we mention, (1) 
transport coefhcients that depend nonanalytically on the 
frequency, in both classical and quantum systems, (2) sin- 
gularities in transport coefficients near continuous phase 
transitions, (3) enhanced, compared to equilibrium, light 
scattering in nonequilibrium systems, and, (4) quantum 
phase transitions that are in different universality classes 
than they would be in the absence of these effects. The 
first three of these effects have been observed. Regard- 
ing the fourth one, the predicted first order nature of 
the quantum phase transition in clean systems has been 
observed, but the critical behavior in disordered systems 
remains to be investigated experimentally. 
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